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ABSTRACT 


The  effects  of  the  transverse  component  of  velocity 
fluctuations  on  nominal  stagnation  point  (NSP)  pressure 
fluctuations  are  analyzed  for  a  finite-thickness  symmetric 
body  penetrating  into  a  turbulent  flow.  The  analysis  is 
idealized  to  the  case  of  a  two-dimensional  symmetric  airfoil¬ 
like  profile.  Considerations  are  limited  to  fluctuation 
scales  which  are  large  relative  to  airfoil  thickness  but  not 
necessarily  large  relative  to  airfoil  chord.  A  single  wave- 
number  component  of  the  turbulence  is  analyzed,  permitting 
synthesis  for  any  desired  wave-number  spectrum,  within  the 
above-mentioned  limitation  on  minimum  scale  size  or  maximum 
wave  number . 

The  results  are  given  in  analytical  form  for  transverse- 
velocity-fluctuation  induced  pressure  fluctuation  at  the  NSP. 

These  results  are  then  compared  to  a  measure  of  axial-velocity 

induced  fluctuations.  The  important  parameters  in  this  com- 

_ .• _ _ _  fluctuation  scale  ,  fluctuation  velocity 

nose  radius  penetration  velocity 

When  the  product  of  these  ratios  isj^  1  20  the  transverse-velocity 

I  is  dominated  by  \  *  \ 


fluctuation  v'  I  is  of  the  same  effectiveness  as  the  axial- 
\  dominates  / 

velocity  fluctuation  u'  in  inducing  pressure  fluctuations  p' 
at  the  NSP.  For  the  case  where  the  axial  component  dominates, 
existing  references  give  the  relation  between  u'  and  p'  adequately. 
For  the  case  where  the  transverse  component  v'  dominates,  the 
present  report  gives  the  relation  between  v'  and  p'.  The 
intermediate  regime  requires  a  combination  of  the  two  contributions, 


T 


I .  •  INTRODUCTION 


This  report  deals  with  certain  aspects  of  the  general 
problem  of  relating  the  pressure  fluctuations  sensed  at  the 
nose  of  a  symmetric  body  to  the  velocity  fluctuations  present 
in  the  fluid  medium  through  which  the  body  is  penetrating. 
Considerations  are  limited  to  low-speed  or  incompressible 
flows.  The  experimental  work  of  Bearman,!  together  with  the 
analytical  background  cited  therein,  forms  a  basis  for  relat¬ 
ing  the  two  aforementioned  fluctuation  quantities  when  the 
scales  of  the  encountered  fluctuations  are  much  smaller  than 
the  body  length  in  the  streamwise  direction. 

Here,  however,  we  are  interested  specifically  in  situa¬ 
tions  for  which  the  encountered  fluctuation  spectrum  (in  wave- 
number  space)  has  significant  energy  in  the  wave-number  regime 
corresponding  to  the  chord-length  of  the  penetrating  body. 

In  the  small-scale  regime  considered  by  Bearman,  the  pressure 
fluctuations  could  be  related  roughly  to  the  axial  component 
(or  streamwise  component)  of  velocity  fluctuations.  It  is 
not  clear  that  this  is  still  true  when  the  encountered  fluctua¬ 
tion  scales  become  as  large  as  the  chord  length  of  the 
penetrating  body,  since  the  transverse  velocity  fluctuations 
then  start  to  affect  the  nose-region  flow  via  the  nonsteady 
circulation  induced  about  the  body. 

Thus  the  objective  of  the  present  report  is  to  estimate 
the  significance  of  the  transverse  component  of  encountered 
velocity  fluctuations  upon  pressure  fluctuations  sensed  at 
the  nominal  stagnation  point  (NSP)  of  a  symmetric  airfoil-like 
profile  traversing  a  turbulent  fluid  medium.  More  specifically, 
the  intent  is  to  determine  when,  if  ever,  the  transverse 
fluctuation  component,  through  the  mechanism  of  inducing 
stagnation-point  meandering,  can  become  as  important  as  the 
axial  fluctuation  component  in  affecting  the  fluctuation 
pressure  sensed  at  the  NSP. 

In  view  of  the  complexity  of  a  completely  detailed 
analysis,  the  approach  to  be  taken  is  as  follows:  A  symmetric 
two-dimensional  wing  moving  at  velocity  UQ  at  zero  nominal 
incidence,  is  presumed  to  penetrate  a  sinusoidal  gust  of 
transverse  velocity  v(x-UQt)  in  wing-fixed  coordinates.  That 
is, the  velocity  fluctuation  is  assumed  to  be  essentially  sta¬ 
tionary  in  fluid  coordinates  or  to  be  convected  in  at  free- 
stream  velocity  U0  in  wing-fixed  coordinates.  This  is  the 
classical  2-D  sinusoidal  gust-entry  problem.  The  convection 
assumption  is  equivalent  to  ignoring  any  mean  flow  in  the 
ambient  fluctuating  field. 


The  present  analysis  goes  beyond  the  classical  gust-entry 
problem  insofar  as  the  classical  problem  treats  a  zero-thick¬ 
ness  flat  plate,  whereas  we  must  not  only  include  thickness 
but,  most  important,  we  must  examine  the  very  nose  point  where 
the  classical  theory  becomes  singular.  Thus  we  must  somehow 
supplement  or  alter  the  classical  gust  entry  problem  to  permit 
etailed  examination  of  nose-region  pressure  fluctuations  for 
an  actual  finite-thickness  profile. 

The  thick-wing  corrections  to  thin-airfoil  theory  that 
are  available  in  the  literature  (to  the  present  author's 
knowledge)  are  limited  to  symmetric  flow  over  symmetric  pro¬ 
files  or  to  steady  asymmetric  flow  (cambered  profiles  or 
Pr°files  at  nonzero  incidence).  6»7  The  present  requirement  is 
for  nonsteady,  non symmetric  flow,  and  a  treatment  will  be 
developed  herein  specifically  for  this  purpose.  The  results 
are  then  compared  with  the  steady,  asymmetric,  thick-wino 
results  of  Reference  7. 

An  approximate  measure  of  the  NSP  fluctuation  pressure 
due  to  the  axial  component  of  velocity  fluctuation  will  be 
provided  as  a  basis  for  comparison  of  the  relative  importance 
of  the  transverse-component  effect,  We  will  be  interested  in 
wavelengths  X_  of  encountered  fluctuations  of  the  order  of  the 
chord  length  Tc)  of  the  penetrating  airfoil.  The  thickness- 
to-chord  ratio,  of  the  airfoil  is  small  of  the  order  of  (e) 
(this  will  be  defined  more  carefully  later) ,  and  it  will  be 
assumed  that  X  is  considerably  larger  than  ec,  i.e.,  that  e 
is  small  and  tnat  the  fluctuation  scales  encountered  are  large 
compared  to  airfoil  thickness.  This  is  a  condition  of  definite 
physical  interest. 

.  ,  1 n  Yiew  of  the  fact  that  we  are  interested  only  in  details 

nose  region,  we  shall  be  content  with  performing 
the  detailed  analysis  for  a  classical  symmett ic  Joukowski 

Pr?f'ie;u  Then'  at  the  end  of  Section  VII,  it  is  shown  how  to 
extend  the  results  to  more  general  symmetric  profiles  by  relat¬ 
ing  the  present  results  to  those  of  Lighthill's  theory  for 
rendering  approximate  solutions  uniformly  valid. 

Intuitiyely ,  since  thin-airfoil  gust-en+ry  theory  leads  to 
a  pressure  singularity  at  the  wing  leading  edge  and  fails 
altogether  to  provide  for  a  stagnation  point,  we  expect  that 
the  effect  of  transverse  velocity  fluctuations  on  stagnation- 
point  oscillation  and  consequent  pressure  fluctuation  at  NSP 
will  become  more  severe  as  the  airfoil  thickness-to-chord 
ratio  (order  e)  decreases.  This  will  turn  out  to  be  true. 

We  also  expect,  intuitively  that  the  effect  will  depend  upon 


4 


the  ratio  of  chord  to  encountered  wavelength  (c/X  ) ,  and  this 
will  be  true  as  well.  What  our  intuition  cannot  predict  is 
the  strength  or  power  to  which  these  effects  enter  nor  the 
appropriate  coefficients,  and  these  will  emerge  directly. 

Summarizing  the  assumptions  to  be  invoked,  we  list  them 
here  concisely  for  convenience  : 

Two-dimensional  approximation 

Symmetric  Joukowski  profile  at  zero  incidence.  It  is 

then  shown,  at  the  end  of  Section  VII 
how  to  generalize  to  more  arbitrary 
symmetric  profiles. 

Zero  mean  velocity  in  "turbulent"  flow  encountered. 

Spatially  sinusoidal  transverse  velocity  fluctuations 

encountered  (these  can,  of  course, 
be  synthesized  into  a  full  spectrum 
of  all  wavelengths,  subject  to  the 
following  restriction) . 

Encountered  wavelengths  \Q  large  compared  to  airfoil 

thickness  but  not  necessarily  large 
compared  to  chord  length. 
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II.  ACCOUNTING  FOR  THICKNESS  BY  MAPPING 

We  take  account  of  airfoil  thickness  by  performing  a 
sequence  of  two  conformal  mappings  from  the  physical  (z) 
plane.  These  are  illustrated  schematically  in  Figure  1. 
First  we  map  from  z,  where  we  have  our  physical  Joukowski 
symmetric  profile,  to  the  intermediate  5-plane  wherein  the 
maps  to  the  (outer)  eccentric  circle.  The  mapping 
is  given  by 


1  £2 

2  =  2  (C  +  C~  (2-1) 

and  the  aforementioned  circle  in  the  5-plane  can  be  given,  in 
terms  of  the  angle  parameter  0,  as 


Cc  =  -  tc  +  t(l+c)elG  (2-2) 

Note  that  angle  0  is  defined  with  respect  to  the  displaced 
center  (-£e,o)  and  not  the  origin  in  the  5  plane.  As 
indicated  on  the  figure,  the  trailing  edge  of  the  physical  air¬ 
foil  in  the  z-plane  is  located  at 


z  =  (£#0) 


(2-3) 


and  the  leading  edge  is  approximately  located  at 


z  =  (--t  (1+2C 2 ) ,  o ) 


(2-4) 


for  small  e. 


The  points  on  the 


physical  profile  corresponding  to 
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Y 


B=Ml+e) 


X 


Z -plane 


Figure  1 


MAPPINGS 


0  =  ±  tt/2  are  given  by 


+ 
z — 


-~l€  +  Uc 


(2-5) 


indicating  a  local  thickness  of  2Jte.  However,  the  maximum 
z-plane  thickness  occurs  (Ref.  1)  for  cosB  =  -  1/2  and  is 
larger  by  a  factor  of  3/T/4.  Thus  the  thickness-to- length 
ratio  is  given  closely  by 


2tc  x 


3/3 

4 


21 


1.3  C 


(2-6) 


Thus  e  has  the  physical  significance 


€ 


~  1 
1.3 


thickness 

chord 


thickness 

chord 


(2-7) 


The  second  mapping,  from  (  to  Z,  is  given  by 


Z  =  (C+U  + 


2  2 
Z+lc 


(2-8) 


and  maps  the  contour  i;,  to  a  slit  of  zero  thickness  in  the 
Z  plane  with  leading  edge  at  (-B,0)  and  trailing  edge  at 
(+B  ,0)  . 

Here 


B  =  Ml+O 


(2-9) 
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is  the  semichord  of  the  Z-plane  slit. 

The  motivation  for  this  double-mapping  is  as  follows: 

We  can  do  the  gust  entry  problem  for  the  slit  (in  the  Z  plane) 
by  classical  techniques.  Then  by  going  back  through  the 
mapping  from  Z  to  5  we  just  cancel  out  the  expected  Z-plane 
leading-edge  singularity.  Then  by  mapping  from  5  to  the 
physical  z-plane,  we  leave  the  singularity  of  the  (5  to  z) 
transformation  inside  the  physical  airfoil  and  retrieve  the 
flow  in  the  physical  plane.  This  sequence  permits  us  to 
utilize  classical  sinusoidal  gust-entry  theory  for  flat 
plates,  and  despite  the  leading-edge  pressure  singularity 
occurring  therein,  transform  the  result  to  a  finite-thickness 
profile  and  compute  the  leading  edge  pressure  fluctuation  on 
this  profile. 

The  requirement  that  the  fluctuation  scale  be  large  with 
respect  to  the  airfoil  thickness  is  necessitated  by  our 
classical  gust-entry  treatment  of  the  Z-plane  flat-plate 
airfoil.  Only  under  this  scale-size  assumption  is  it  legal  to 
assume  that  the  encountered  transverse  velocity  is  the  same 
on  upper  and  lower  surfaces  of  the  physical  airfoil  at  the 
same  chordwise  point,  x.  This,  in  turn,  is  implicit  in  the 
classical  treatment  of  the  Z-plane  slit  where  the  transverse 
gust  to  be  counteracted  by  a  potential  flow  is  identical  on 
upper  and  lower  surfaces  of  the  slit  at  the  same  value  of  X. 


III.  FLOWS  IN  THE  VARIOUS  PLANES 

Defining  the  complex  potential  function  by  F (z)=F(C)=F (Z) 
at  corresponding  points ,  where 


F  “  <p+i  0 


(3-1) 


and  where  4>  is  the  velocity  potential  and  1 1>  is  the  stream  func 
tion,  we  have  for  the  basic  (zero-disturbance)  flow  in  the 
C-plane , 


Po(0  -  f  <C+«  +  )  (3-2) 


In  the  Z  plane,  this  is  simply 


F  (Z)  =  U  Z  (3-3) 

o  o 

which  is  the  purely  uniform  flow  in  the  positive  X  direction, 
as  it  should  be. 

Inverting  Eq.  (2-1)  for  5,  we  find 

C  =  z  +  J  z2-l2  (3-4) 


This  gives,  for  the  basic  flow  in  the  z  plane, 


Po(z) 


u  / - 

“  (z+t€  +  JZ2-t2  + 


l2(l+€)2 


Z+<t  €  +v z  ^  ^ 


(3-5) 


In  the  far  field  (|z|  large)  this  approaches 


Fo(z)  ~  Uc  z 


(uniform  flow) 


(3-6) 


as  it  should.  FQ(z)  is  the  complex  potential  for  uniform  flow 
of  velocity  UD  in  the  positive  x  direction  over  the  Joukowski 
profile  of  Figure  1,  i.e.,  in  the  physical  plane. 

On  the  profile  5=5C  in  the  5-plane, 


Ud+Oe16*-^)' 


l  (1+c )  e‘ 


=  UQ  l  (l+o )  cos0  =  pure  real 


(3-7) 


so  that  ifi=o,  and  the  profile  is  the  zero-streamline  in  all 
three  planes. 

The  complex  velocity 

w  =  u  -  iv  is  given,  in  each  plane  by 

,  .  dF(z) 
w(z)  =  ±~ 

CLZ 

«<«>  =  i§ 

w(z)  =  fl 

(3-8) 


We  can  therefore  find  the  z  (physical) -plane  flow  from 
the  Z-plane  flow  using  the  relation 


li 

.  0 

Id 

ID 


w(z)  -  w(Z) 


=  w  (Z) 


1  - 


iinur 

(£K02 


.  il 
c2 


(3-9) 


At  the  leading  edge  (corresponding  toC=-*(l42E)  the 
denominator  dz/dC  is  regular,  nonzero,  but  the  numerator 
dZ/d5  vanishes.  However,  it  will  turn  out  that  our  required 
nonsteady  flow  in  the  Z  plane  will  approach  infinity  at  the 
leading  edge  (Z-(-B,0))  in  just  such  a  way  as  to  balance  the 
vanishing  of  dZ/d£,  and  thereby  contribute  the  bounded, 
nonzero  velocity  at  the  NSP  in  the  physical  plane. 

Now  in  the  inflow  in  the  Z-plane  we  superimpose  on  the 
uniform  U  -flow  a  frozen  sinusoidal  transverse  gust  whose 
form,  on  the  ^cixis*  is 


with 


v  =  vq  cos(wQt  -  kQX) 

v  <  c  U 
o  o 


CO 

*  =  ^  =  -2 

0  Xo  Uo 


(3-10) 


(3-11) 


■f  00 


together  and, for  large  negative  X, 


Since,  |z|  ,  — |z|  all 


x  ~  X 


(3-12) 


(see  Eqs.  (2-1)  and  (2-8), 

and  since  the  mappings  are  conformal  (angle -preserving),  there¬ 
fore  this  sinusoidal  transverse  gust  must  correspond  or  map  to 
the  same  transverse  inflow  in  the  field  far  upstream  of  the 
wing  in  the  physical  z-plane.  We  must  then  find  that  potential 
flow  which  exactly  counteracts  the  gust  (  Eq.  (3-10))  on  the 
slit  airfoil  in  the  Z-plane.  This  flow  when  mapped  to  the 
z-plane,  will  then  leave  the  physical  airfoil  a  streamline  and 
hence  provide  the  proper  flow  in  the  physical  plane,  when  the 
(far  upstream)  inflow  in  the  physical  plane  is  as  required; 
namely  the  uniform  flow  UQ  plus  the  frozen  transverse  sinusoidal 
gust.  The  solution  for  the  counteracting  potential  flow  in  the 
Z  plane  is  the  classical  thin-airfoil  sinusoidal  oust  entry 
problem,  except  that  we  need  the  vorticity  or  pressure- jump 
distribution  in  detail,  whereas  the  usual  derivations  are  content 
with  integrals  for  lift  and  moment. 


IV. 


PRESSURE  AT  NSP  IN  THE  PHYSICAL  PLANE 


After  we  have  found  the  required  counteracting  potential 
flow  we  must  then  utilize  it  to  construct  the  pressure  at  the 
nose  (NSP)  of  the  physical  airfoil.  For  this  we  use  the 
nonsteady  Bernoulli  equation 


[if  +  p  +  poqJ/2].  -  p  +  °o^/2  *  po  f?  (4-1) 

Now,  as  noted  (Eq.  (3-10))  in  the  preceding  section. 


v  «  U  or  v  <  cu 

o  o  o  o 


so  that  the  left  side  of  (4-1)  is  essentially  equal  to 


p.*P  U  2/2  + 

00  o  o  at 


+  0(Uo2C2po) 


(4-2) 


On  the  richt-hand  side  of  Eq.  (4-1)  we  need  the  velocity 
q  **  |u+iv|  which  consist*.'  of  superimposed  contributions  from 
the  undisturbed- flow,  the  sinusoidal  inflow,  and  the  gust¬ 
counteracting  potential  flow.  At  the  NSP,  the  first  of  these 
vanishes,  since  NSP  is  a  stagnation  point  of  the  undisturbed 
flow  (Eq.  (3-5)).  The  second  contribution  (the  inflow  sinu¬ 
soidal  gust)  vanishes  also  since  it  is  finite  (small,  actually) 
in  the  Z-plane  (Eq.  3-10)  and  is  annihilated  by  dZ/d t  at  the 
leading  edge  or  NSP,  in  accordance  with  Eq.  (3-9)  and  the 
remarks  following  that  equation. 

without  closer  investigation,  it  would  appear  that  the 
third  contribution  to  q  at  the  NSP  would  vanish  for  the  same 
reason.  However,  as  indicated  earlier,  this  component  (the 
counteracting  flow)  proves  to  be  singular  of  just  the  correct 
strength  (as  X  -*•- B)  as  to  offset  the  vanishing  of  dZ/dC  and 
contribute  to  q  at  the  physical  NSP.  We  shall  devote  most  of 
the  remainder  of  the  report  to  this  computation. 


The  terms  3$/3t  in  Eq.  (4-1)  should  be  examined  (as 
they  may)  in  the  Z  plane.  The  contributions  to  3$/3t  on 
the  right-hand  side  of  Hq.  (4-1)  come  from  the  inflow  gust  and 
the  counteracting  nonsteady  flow.  But  the  first  of  these  is 
exactly  balanced  by  the  3$/3t  on  the  left-hand  side.  We  are 
assuming  here  that  the  sinusoidal  gust  inflow  can  be  described 
by  a  potential  flow.  This  is  not  true  of  physical  turbulance 
but  it  is  sufficient  for  the  present  single-sinusoidal  argu¬ 
ment,  and  we  can  easily  construct  potential  flows  having  the 
prescribed  transverse  form  (3-10)  on  the  X  axis. 


Thus  all  that  survives  of  the  3$/3t  terms  on  both  sides 
of  Eq.  (4-1)  is  the  contribution  to  the  right  hand  side  due 
to  the  gust-counteracting  flow.  It  will  turn  out  that  this 
vanishes  on  the  upstream  X-axis  right  up  to  the  leading  edge 
and  hence  vanishes  at  the  NSP  in  the  physical  plane.  This  will 
emerge  from  the  gust-entry  calculation  but  is  actually  a 
direct  consequence  of  the  fact  that  the  Z-plane  gust-entry 

iS  odd.in  Y'  vanishinq  unstream  of  the  silt  and 
exhibiting  an  antisymmetric  jump  on  and  downstream  of  the  slit. 

Since  the  jump  in  $  is  regular  (nonsingular  at  the  leading 
edge)  therefore  $  and,  with  it  3^/3 t, vanishes  at  the  NSP. 

The  consequence  of  the  foregoing  statements  with  respect 
to  the  pressure  fluctuation  at  the  NSP  is  as  follows: 


P  - 


(4-3) 


wherein  q.  is  the  velocity  magnitude  due  to  the  gust-entry- 
counteracting  potential  flow  in  the  Z  plane,  mapped  to  the 
NSP  of  the  z  plane.  That  is 


^1 


C  U(z)-iv(Z)  ] 
Z-f-B.O1) 


dZ 

dC 

dz 

dC 


(4-4) 


V.  LIMITING  BEHAVIOR  OF  TELOCITY 
TRANSFORMATION  AT  NSP 

We  examine  in  detail,  in  this  section,  the  limiting 
behavior  of  the  velocity  transformation  equation  (3-9)  as 
the  NSP  is  approached  in  the  physical  plane  or,  what  is 
equivalent. 


as 


or 


Z  -  (-B  ,0) 


z  -  (-£(1+20  ,0)  / 

Inverting  Eq.  (2-8)  we  find 

\ 


+  Z  +  / z2-b: 


C  “  “ 

where  B  =  £(1+0 


) 


(5-1) 


(5-2) 


We  consider  points  Z±  on  the  top  and  bottom,  respectively,  of 
the  slit  and  accomplish  the  approach  to  the  leading  edge  as 
follows:  We  let 


Z^ 


B  +  r 


(5-3) 


top 

bottom 


with  r  small  with  respect  to  cl  and  positive,  real. 


r  «  c  £ 


(5-4) 


Then,  from  (5-2) 


corresponding  points  £*  are  given  by 


£c  -  B  +  r  + 


i 


/ 


2Br-r2 


(5-5) 


14 


II 


D 

I! 

li 

li 

I! 

II 

II 

I! 

0 

II 


Then  the  denominator  dz/dc  =  (1-  -7-) in  Eq.  (3-9)  remains 
finite  as  r*o.  ^ 


dz 

dC 


r~o 


1  - 


2  5f 

*T(l+2c) 


1- 


1 

(l+2e) 2 


=  i  -  [l-4c-4c2+16c2 ] 


*  4C  (1-3C) 


(5-6) 


Since  our  Z-plane  gust-entry  analysis  will  be  good  only 
to  order  e,  therefore  we  keep  only  terms  to  this  order  in 
Eq.  (5-6)  and  write 


dz 

dC 


r-o 


4c 


(5-7) 
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The  numerator  dZ/dG  in  Eq.  (3-9)  is  more  involved.  It 
can  be  written  as 
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Keeping  only  terms  of  order  and  utilizina 

results  (5-6)  and  (5-9)  in  Eq.  (3-9)  we  find 


vr(z) 


Z=(-B+r) 


top 

bottom 


(5-10) 


Now  w  =  u-iv,  and  the  counteracting  part  of  v  simply  cancels 
out  the  sinusoidal  gust  of  Eq.  (3-10),  which  is  finite  every¬ 
where.  Thus  the  v  contribution  to  (5-10)  is  annihilated  as 
r-*o,  i.e.  as  we  approach  the  nose  or  NSP  in  the  physical  plane. 

However,  the  u-part  of  w  on  the  slit  is  singular.  In 
fact,  foreseeing  the  form  in  which  the  result  will  emerge,  we 

WT*  i  *  “ 


u  (X  ,  0~)  = 


B+X 


(5-11) 


where,  in  f  act,  f  is  regular  for  all  X  and  depends  harmonically 
on  time  t.  We  simply  Veep  t  as  a  parameter  for  the  present. 

We  will  find  the  function  f  in  the  next  section.  However,  we 
see  that  as 
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X  -  (-B+r) 


u(X,0±)  -  ±  \  f  (-B,  t ) 

Then,  from  Eq.  (5-10' 


or 


u  -  iv 


i 

2C 


f (-B,t) 


(5-12) 


(5-13) 


Since,  physically, f  must  be  real,  we  have 


u  -  o 

q  -  v  -  1L- B,t) 
2c 


as  the  NSP  is  approached  (as  r+o)  . 

This  is  physically  reasonable,  since  the  vertical  tangent  at  the 
nose  (in  the  physical  plane)  prohibits  any  axial  velocity  u  and 
permits  only  the  vertical  component  v.  We  see,  from  Eq.  (5-14) 
how  the  velocity  induced  at  the  NSP  varies  inversely  with  c  ' 
which,  we  recall,  is  0.77  of  the  thickness-to-chord  ratio. 

% 

The  following  section  is  devoted  to  the  development,  in 
detail,  of  the  function  f(X,t)  first  introduced  in  Eq.  (5-11). 

In  this  development,  the  function  sought  will  be  the  time- 
dependent  vorticity  density  function  y(X,t)  which  exists  over 
the  Z-plane  slit  and  its  wake  (from  +b  to  «)  on  the  X  axis. 
(5-11)  is  equivalent  to  foreseeing  the  form 
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y  (x,t) 


f (x,t) 


(5-15) 


where  vorticity  and  circulation  are  designated  positive  for  a 
clockwise  sense  of  rotation 

o 

and  where  f  is  real  and  everywhere  regular. 

It  can  be  shown  that  the  result  of  our  double  mapping 
procedure  (Eq.  (5-10))  for  the  transformation  of  the  complex 
velocity  in  the  nose  region  of  the  slit  (Z)  plane  to  that  in 
the  nose  region  of  the  physical  (z)  plane,  agrees  in  magnitude 
with  the  correction  factor  of  Lighthill's  nonlinear  theory  for 
steady  flow  over  lifting  (asymmetric)  airfoils.  This  follows 
from  the  fact,  as  pointed  out  by  Van  Dyke,'6*  that  the  leading 
edge  radius  (a)  of  a  Joukowski  symmetric  airfoil  is  given  by 


a 


x  semichord 


=  2c2  x  chord 


(5-16) 


Then,  as  noted  by  Thwaites,  the  first-order  Lighthill  correct 
tion  for  a  two-dimensional  profile  of  nose-radius  (a)  is 
achieved  by  multiplication  of  the  linear-theoretical  velocity 
by  the  factor 


(5-17) 


where  s  is  distance  aft  of  the  leading  edge.  In  our  notation, 
with 
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and  this  agrees,  in  magnitude,  with  Lighthill’s  result, 

Eq.  (5-20).  In  Section  VII  herein,  we  use  this  agreement  to 
generalize  our  results  to  shapes  other  than  the  Joukowski 
profile. 
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VI.  .  THE  GUST  ENTRY  PROBLEM  IN  THE  Z-PLANE 
(DETAILS  OF  f'X,t))  - 

Most  of  the  analysis  needed  for  the  gust— entry  problem  for 
a  flat  plate  is  available  in  sources  such  as  Reference  2. 

However,  we  depart  slightly  in  some  areas  from  the  classical 
derivation  to  keep  certain  details  clear  and  to  focus  attention 
on  the  results  required  for  the  present  application. 

The  gust-counteracting  flow  in  the  Z  plane  is  completely 
aescribable  in  terms  of  a  time-dependent  bound  vorticity 
Ya(X,t)  and  a  trailing  or  shed  vorticity  yw(X,t) ,  arising  as  a 

consequence  of  the  time  c  cpenucnce  of  the  total  bound  circulation 
*  v  tz. )  • 


B 

r(t)  =  J*  _  (x,t)  dx 

“13  ci 


(6-1) 


Temporarily  we  admit  complex  time  dependence  ela)°  (which 
leads  to  complex  ya,  yw)  with  the  understanding  that  we  will 
return  to  real  quantities  before  involving  the  mapping  from  Z 
back  to  z .  The  inflow  gust  to  be  counteracted  is  qiven  bv 
(Eq.  (3-10))  1 


v 


i  (co  t-k  X) 
o  o 


(6-2) 


Then  our  counteracting  flow  must  be  given 
tions  y  /  y  such  that  (dropping  the  elut 
is  common  to  v,  Yu/  yw/  and  D 


by  vorticity  distribu¬ 
ting  dependence  which 
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T  B  y  (x')dx' 
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2n  j-b  x'-x 


i_  -  rw(x’)dx- 
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X'-X 


(6-3) 
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The  assumption  (3-10)  that  v0<eU0  is  what  permits  us  to 
approximate  the  slightly  wavy  vortex  wake  geometry  with  the 
X-axis  proper.  It  is  one  of  the  reasons  why  our  gust-entry 
analysis  is  valid  only  to  order  e,  as  was  noted  in  Section  V. 

Scaling  lengths  with  respect  to  B  and  using  (temporarily) 


x_  ■ 


(6-4) 

and  using  the  classical  reduced  frequency  (corresponding  to 
the  slit  in  the  Z-plane) 

co  B 

*  u  (6-5) 


we  rewrite  Eq.  (6-3)  in  the  scaled  form 


+  277  Jl  £-x 


€-x 


(6-6) 


where  ^  denotes  y  in  scaled  variables  ,x,£  . 

Now  the  requirement  of  zero  pressure  jump  across  the 
wake  demands  that,  in  physical  variables. 


y„  ■  (6'7) 


This,  in  turn,  requires  that 
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y  =  e  l  of 

9  T.f 


y  (b) 

w 


or,  in  terms  of  scaled  variables, 


~  _  ik  -ikx  .  . 

yw  -  e  e  >W(B) 


Then,  vorticity  conservation  demands  that 

.  ST  *  (B)  dx 

dt  ywlB)  dt 


-  iw  T  =  U  y  (B) 
o  o  'w 


-  ““  -  y  (b) 

B  w 


(6-10) 


Thus 


YU)  =  eik  e'ikx 

W  B 


(6-11) 


Reference  (2)  defines  a  convenient  parameter  ft: 


B  e 


(6-12) 


so  that 


y  (x)  =  -  ik  n  e 

w 


-ikx 


(6-13) 


The  value  of  ft  will  be  determined  subsequently.  Then  Eq.  (6-6) 
can  be  written 
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2ir  J  £-x 
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(6-14) 


The  inversion  of  this  classical  singular  integral 
equation  for  ya(Z)  is  nonunique  unless  we  invoke  some 
additional  requirement. The  Kutta  condition  is  the  appropriate 
one  and,  in  the  present  context,  takes  the  form  of  a  boundedness 
requirement  on  YaU)  at  the  trailing  edge  (£=1) .  With  this  con¬ 
dition,  the  inversion  is  rendered  unique  and  takes  the  form  (4) 


ya(x) 
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_  ,1 

=  _  1_  /l-x  /  /T+l  1  I 

JT  y  1+x  J  J  l-£  £-x  ]  Vc 
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,-ik« 


.  ilcQ  f  e~ik?7c 

2*  J  H-C 


(6-15) 


Because  n  varies  over  (l,00  )while  £  varies  over  (-1,1),  it  is 
legal  to  invert  the  order  of  the  £  and  n  integrations  (if  they 
were  in  the  same  range  this  lyould  not  be  legal  but  would  re¬ 
quire  the  use  of  the  Poincare -Bertrand  relation  of  singular 
integral  equation  theory).^)  Performing  the  inversion  of 
order,  we  ^ind 


We  introduce  the  substitutions 


£  *=  cos© 
x  *  cos <p 

tj  -  cosh  s 
Then  Eq.  (6-16)  takes  the  form 


(6-17) 
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cosQ-cos  <t>  I  l  cos©-cosh  s 


(6-18) 


The  first  term  gives 


(1) 


IT 

T)-X 

while  the  second  may  be  written  as 


(6-19) 


1  nff  1+  cos9  -  cosh  s  +  cosh  s 
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cos  6  -  cosh  s 


-  <  n  +  (1+cosh  s) 

*/“x  I  o 


cos 0-cosh  s 


Iff-  (1+cosh  s)  — ~r — 
T)-x  (  smh  s 


as  may  be  easily  shown  from  integration  by  residues. 
Thus,  we  find 


+  _^L_ 
Vt?2-1  ; 


r?-x 


=  _JL_  La ±1 . 

7?-x  I  rj-1 


Inserting  this  into  Eq.  (6-15)  we  find 
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Now,  by  definition 
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Thus  we  find  ft  by  operating  on  both  sides  of  Eq.  (6-22)  with 
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)  dx 


(6-24) 


This  gives,  upon  inverting  integration  orders,  as  we  may. 
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dx 
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The  first  inner  integral  gives 
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(6-25) 
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while  the  second  inner  integral  gives 
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Thus 
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The  second  integral  is,  according  to  Ref.  (2)  f 
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giving  the  result 


(k)+i 


H^2) (k) ] 
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ik-i) 


(-voo 


-ikt 


or 
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(6-30) 


Most  of  the  relations  developed  in  this  section  thus  far 
ccn  be  found  directly  in  or  inferred  from  the  flat-plate  gust 
entry  theory  presented  in  Reference  (2) .  They  are  presented 
here  primarily  to  make  the  present  report  self-contained. 

Comparing  the  form  of  Eq.  (5-15)  for  y(X,t)  with  the 
results  of  this  section  (Eq.  (6-22)  together  with  Eq.  (6-30) 
for  fi)  we  see  that  the  regular  function  f(X,t)  is  given  on  the 
wing  (-B ,B)  by 


(6-31) 


with  «  given  by  (6-30) . 
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But  all  we  need  for  the  present  application  (Eq.  (5-14)) 
is  f(-B,t).  For  this  we  can  use  (6-31)  with  X  set  equal  to 

-B. 

Thus 


(6-32) 
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The  first  integral  in  (6-32)  gives 
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The  second  integral  in  (6-32)  gives 
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(6-34) 


(6-35) 


The  integral  required  in  the  expression  for  ft  (Eq.  (€-30)) 
may  be  performed,  with  the  result,'3* 
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Finally,  we  find  for  f(-B,t),  from  (6-35) 

it*  \  kH(2)(k)  .. 

f(-B,t)  =  eL  ±  1  2  v  J  (k)  -  — 2- -  •  ^  v  * 

1  o  o'  2  k  o 


JQ  (k ) - i (k) 

H*2*  (k)+iH*2*  (k) 
l  o  * 


33 


ia*  J  2(H^2)+1H^2)-h'2))(J  -iJ  ) 

6  Vo  2  Jc<k>  -  ~  \a,  °  (2 ) 1 - -0— 

(H^[  '  +  i  IT  ') 


ioot  J  .  ^  Hj  ^k)(Jo(k)-iJl(k;)) 

e  2  v  /  l  J  (k)  +  j2)  (2) 

Bpk)  +  iH^'(k) 


Now  the  quantity 
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is  the  famous  complex 


Theodorsen  C(k)  function  of  oscillatory  thin-wing  theory. 
Nov; 

co  B  co  c  U  k  c 
o  ~  o  o  o 


c  2t r  rrc 


(6-39) 


since 


B  =  f 


Thus,  for  inflow  gust  forms  having  wavelength  Xc  of  the  order  of 
chord  c,  we  are  in  the  region  where  k  is  of  the  order  of  tt,  which 
is  large  in  oscillatory  aerodynamic  theory.  For  such  large  k. 


C(k)  =  0.5  +  i  x  0=0.5 


and  f(-B,t)  simplifies  to 
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If  we  make  the  further  asymptotic  approximation  for  the 
Jo00  /  J^(k)  as  well,  (this  is  not  as  well  justified  with  k  of 
the  order  of  it,  but  if  anything  it  slightly  underestimates  the 
contribution  considered)  then 
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Returning  to  real  form  and  unsealed  variables. 


(5-42) 


f(-B,t)  = 


=  vo  yVirr  cos<wt  +  k0  B  -  !> 


(6-43) 


Inserting  this  into  our  relation  (5-14)  for  q,  at  the 
NSP  in  the  physical  plane  1 
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For  the  case  where  encountered  fluctuating  scales  are 
markedly  larger  .than  the  chord  length  of  the  penetrating  airfoil 
then  the  simple  asymptotic-approximation  form  (Eq.  (6-44))  for 
should  be  replaced  by  the  more  accurate  relation  arising 
from  Eqs.  (5-14)  and  (6-38);  namely, 
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The  approximation 
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Hl(2,+  i  H 
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=  0.5  +  i  (o) 


is  quite  good  for  k  values  as  small  as  1.0. 

The  asymptotic  approximations  toJQ  and  Ji,  however,  should 
be  replaced  by  their  actual  respective  values  for  situations  in 
which  encountered  fluctuation  scales  far  exceed  the  penetrating 
airfoil  chord  length  and  k  is,  accordingly,  less  than'-Ti. 

In  support  of  our  statement,  in  Section  IV,  concerning 
the  vanishing  of  S^/St  at  the  NSP,  we  note  that  the  complex 
potential  function  for  our  gust-counteracting  flow  is  given  by 


F(Z)  =  ifl  +  i  |)  = 
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For  real  (physical)  y(X,t)  this  gives 


y(x',t)  tan-1  dx' 
x-x' 


(6-47) 


and  we  see  that  4>  and  3<f>/3t  are  odd  in  Y  and  therefore  vanish  on  the 
X  axis  upstream  of  X=-B.  Since  the  jump  in  <J>  across  the  airfoil  is 
not  singular,  in  contrast  to  y(X,t),  therefore  $  and  3<j>/3t  vanish 
at  the  nose  (X=-B)  and  therefore  vanish  at  the  NSP  in  the  physical 
plane  as  well. 

Thus  we  have  finally  for  the  fluctuating  pressure  at  the 
NSP  in  the  physical  plane 


v2  cos2  (cot+k  B-“) 
o  o  4 


2C  It  k  B 
o 


(6-48) 
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VII 


.  COMPARISON  OF  PRESSURE  FLUCTUATIONS 

INDUCED  BY  TRANSVERSE  AND  LONGITUDINAL 
VELOCITY  FLUCTUATIONS 

As  an  approximate  measure  of  the  effect  of  pressure 
fluctuations  induced  by  axial  velocity  fluctuations 


u  =  U  +  6u  cos  (o>  t-k  x) 
o  o  o 

we  simply  take  the ‘quasi-steady  result  * 


(7-1) 


fip 


U 

o 


fiu 


*due  to  6u 


(7-2) 


Assuming  6u  of  the  same  order  as  v  (isotropy  in  planes  normal 
to  the  wing  span)  we  see  that  the  ratio  y  of  pressure  fluctuation 
amplitudes  induced  by  transverse  velocity  fluctuations  to  those 
induced  by  axial  fluctuations  is  given  approximately  by 


v 

o 


ff  U  k  B 
o  o 


v  X 
o  o 


2  2 

8  € it  B  U 


(7-3) 


Recalling  that  e  = 


thickness 

chord 


f 


chord 


.V 


1  thickness  J 


Since  B  = 


2  chord , 


this  is  equivalent  to 


(7-4) 


_  (5) 

*In  fact,  Bearman's  experimental  results  with  a  bluff 
body  show  this  to  be  satisfied  within  5%  for  scales  XQ  more 
than  five  times  as  big  as  body  thickness,  and  we  are  restrict¬ 
ing  our  treatment  to  this  regime. 
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Vp  _ chord  27 

U  thickness  thickness  64 
o 


chord 


1 


thickness  thickness ^ 


(7-5) 


This  makes  clear  what  is  required  to  give  physical  sig¬ 
nificance  to  the  transverse-gust-induced  contribution  to  pressure 

fluctuations;  namely  large  -  — ■ -g — —  ,  large 

thickness  thickness 

and  slow  penetration  velocity  UD  or  strong  fluctuation  intensity 
(vQ)  . 


If  Xo  ~  chord 

j  chord 
ana  . .  >  p—  ~~ 
thickness 

Then  M  = 


~  10 
100  Xo 

23  U 


(7-6) 


while  for 


V 

_n  ~ 
uo 

1% 

~  L_ 

23 

vo 

-  -  -- 

10% 

,  u 

~  _L. 

Uo 

2.3 

(7-7) 


We  should  now  be  able  to  remove  any  restriction  implicit 
in  our  use  of  a  Joukowski  profile  by  replacing  the  ez ,  which 
appears  in  expression  (7-3)  for  yf  by  the  relation  (Eq.  (5-16)) 
between  e2  and  nose-radius  (a) 


c2  =  a _  (7-8) 

2  x  chord 

(7) 

Then,  in  accord  with  the  Lighthill  theory,  we  can  write  our 
ratio  y  in  terms  of  profile  nose  radius  and  thereby  add  generality 
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to  the  result.  Thus,  in  terms  of  nose  radius,  for  any  symmetric 
profile  of  small  thickness-to-chord  ratio. 


v  X 
o  o 


chord 


B  U 


v  X 
o  o 

7T 


2a  irz  U 


since  B  = 


c 

2 


20  U  nose  radius 
o 


(7-9) 


This  expresses  the  relative  importance  of  transverse  and 
axial  fluctuation  velocity  components,  in  affecting  NSP  pressure 
fluctuations,  directly  in  terms  of  the  ratios  of 
fluctuation  velocity  fluctuation  scale 

penetration  velocity  ancJ  nose  radius  * 

For  a  profile  whose  nose  region  resembles  that  of  a  symmetric 
Joukowski  profile,  expressions  (7-9)  and  (7-5)  should  agree.  If 
the  nose  region  of  the  profile  of  interest  differs  markedly  from 
the  Joukowski  nose  profile,  then  relation  (7-9)  should  still  be 
valid. 


VIII. 


CONCLUSIONS 


The  results  of  our  analysis  relating  NSP  pressure  fluctua¬ 
tions  to  encountered  velocity  fluctuations  are  contained  in 
expressions  (6-48)  ,  (7-5)  ,  and  (7-9) .  We  reiterate  that  the 
present  analysis  applies  only  to  fluctuation  scales  which  are 
large  relative  to  profile  thickness.  They  may,  however,  be  of 
the  same  order  as  profile  chord-length.  The  aforementioned 
results  indicate  that  as  the  ratios 

fluctuation  scale  _  ^o 
nose  radius  “  a 


and 

fluctuation  velocity  _  ^o 

penetration  velocity  UQ 

increase  such  that  their  product  approaches  a  number  of  the 
order  of  20,  then  transverse  velocity  fluctuations  become  as 
important  as  axial  fluctuations.  For  fluctuation  scales  whose 
length  is  at  least  an  order  of  magnitude  less  than  the  profile 
chord,  Bearman ' s  results should  serve  to  relate  NSP  pressure 
fluctuations  to  encountered  axial  velocity  fluctuations.  For 
fluctuation  scales  of  the  order  of  profile  chord  or  longer, 
Bearman' s  results  probably  still  apply  if  conditions  are  such  that 
our  transverse  velocity  effectiveness  ratio  y  is  still  small. 

This  would  require  verification.  When  y  is  not  small,  then 
Bearman' s  results  cannot  be  used,  since  they  neglect  the  effects 
of  transverse  velocity  altogether. 

When  conditions  are  such  that  y  is  large,  then  the  trans¬ 
verse  velocity  fluctuations  dominate  the  fluctuation  pressure  at 
the  NSP,  and  our  expression  (6-48)  gives  the  result.  As  in 
Section  VII  herein,  we  can  add  generality  to  Eq.  (6-48)  (which  was 
derived  for  a  Joukowski  profile)  by  replacing  e  therein  with 


r  a  — 

V  2  x  chord 


With  this  substitution  our  expression  (6-48)  should  give  the 
pressure  at  the  NSP,  valid  for  the  large-y  regime,  for  nose 
shapes  more  general  than  that  characterizing  the  Joukowski  profile. 

For  the  case  where  y  is  large  and  the  encountered  fluct¬ 
uation  scales  are  large  with  respect  to  penetrating  airfoil  chord 
length,  then  the  approximate  Eq.  (6-44)  for  q^  should  be  replaced 
by  the  more  accurate  expression  (6-45)  in  forming  the  fluctuation 
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pressure  at  the  NSP.  Equation  (8-1),  replacing  e  by  its 
equivalent  in  terms  of  nose-radius,  can  still  be  used  within 
this  improved  relation  (6-45)  to  remove  any  restrictions 
imposed  by  the  assumption  of  a  Joukowski  profile. 
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